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Abstract 

A permutomino of size rt is a polyomino determined by a pair 
(''■i,7r2) of permutations of size n + 1, such that TTi{i) ^ 7i'2(*), for 
l<i<n+l. In this paper, after recaUing some enumerative results 
about permutominoes, we give a first algorithm for the exhaustive gen- 
eration of a particular class of permutominoes, the convex permutomi- 
noes, proving that its cost is proportional to the number of generated 
objects. 



1 Introduction 

A permutomino is a special polyomino, defined by two permutation matrices 
having the same size. The class of permutominoes was introduced by Incitti 
in [To] while studying the problem of determining the i?-polynomials (re- 
lated to the Kazhdan-Lusztig R-polynomials) associated with a pair (/i, v) 
of permutations. In his paper Incitti gave a general definition of these com- 
binatorial objects using some algebraic notions. 

In this paper we use the definition of permutominoes given in [7] which 
does not use the algebraic notions and nevertheless, though different, it turns 
out to be equivalent to Incitti's one. 

The main results about permutominoes concern the enumeration of var- 
ious subclasses of permutominoes and the characterization for the permuta- 
tions defining these subclasses [31 EJ El EJ [9] , while at our knowledge nothing 
exists about their generation. On the other hand, exhaustive generation of 
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combinatorial objects [Tl[31[5] is an area of increasing interest. In fact, many 
practical questions in diverse areas, such as hardware and software testing, 
and combinatorial chemistry, require for their solution the exhaustive search 
through all objects in the class. 

Actually, in [7] a recursive generation of all convex permutominoes of 
size (n + 1) from the ones of size n, according to the ECO method [2], 
is presented. Section [2] contains basic definitions and some enumerative 
results of convex permutominoes of size n. Section [3] recalls the recursive 
generation of convex permutominoes presented in [7] and Section[l]illustrates 
the exhaustive generating algorithm based on the recursive construction 
recalled in Section [3l 

2 Basics on permutominoes 
2.1 Definitions and properties 

In order to define permutominoes we need to introduce polyominoes. In the 
plane Z x Z a cell is a unit square and a polyomino is a finite connected union 
of cells having no cut point. Polyominoes are defined up to translations. A 
column {row) of a polyomino is the intersection between the polyomino and 
an infinite strip of cells whose centers lie on a vertical (horizontal) line. 




(a) (b) (c) 



Figure 1: (a) a column convex polyomino; (6) a convex polyomino; (c) a 
directed (not convex) polyomino with a hole. 

In order to simplify many problems which are still open on the class 
of polyominoes, several subclasses were defined by combining two notions: 
the geometrical notion of convexity and the notion of directed growth. A 
polyomino is said to be column convex [row convex] if its intersection with 
any vertical [horizontal] line is convex (Figured] (a)). A polyomino is convex 
if it is both column and row convex (Figured] (6)). In a convex polyomino 
the semi-perimeter is given by the sum of the number of rows and columns, 
while the area is the number of its cells. 
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A polyomino P is said to be directed when every cell of P can be reached 
from a distinguished cell, called root (usually the bottom leftmost cell), by 
a path which is contained in P and uses only north and east unit steps 
(Figure [11(c)). 

Let P be a polyomino without holes having n rows and n columns, n > 1; 
without loss in generality, we assume that the bottom leftmost vertex of the 
polyomino minimal bounding rectangle lies in (1,1). Let ^i, . . . , ^42(^+1) 
be the sequence of the vertices of P obtained by visiting the boundary in 
clockwise sense, starting from its leftmost point with minimal ordinate. 

We say that P is a permutomino if the sets Vi = {^1,743, . . . ,A2r+i} 
and V2 = {^2; ^4; ■ ■ ■ ) ^2r+2} represent two permutations matrices of [n + 1] 
= {1,2,. . . , n + 1}. Obviously, if P is a permutomino, then r = n, and n 
is called the size of the permutomino. The two permutations defined by Vi 
and V2 are indicated by (7ri(P), 7r2{P) ) (briefly, (7ri,7r2)), respectively (see 
Figure [2]). Given a pair of permutations to = (/x, of [n + 1], we say that 
a permutomino P is associated with if /i = 7ri(P) and v = ■K2{P)- 



71,^(1, 3, 4, 2, 5) 
It ,= fi, 4, 5. J. 2) 



71 , = fi, 3, 4, 2, 5) 
n2=(3, 4, 2. 5. 1) 



Figure 2: Two permutominoes and the associated permutations. The per- 
mutation TTi (resp. 712) is represented by black (resp. white) dots. 



A permutomino is convex {directed) if it is a convex (directed) polyomino. 
A parallelogram permutomino is a directed and convex one having the (1, 1) 
and (n, n) vertices in common with its minimal bounding square; a stack 
permutomino is a directed and convex one in which the bottom side of its 
minimal bounding square belongs to the permutomino itself. 

The definition of permutominoes leads to the following remarkable prop- 
erty: 

Proposition 1 Any permutomino P has the property that, for each abscissa 
(ordinate) there exists exactly one vertical (horizontal) side in the boundary 
of P having such coordinate. This property is also a sufficient condition for 
a polyomino to be a permutomino. 

Starting from the leftmost point having minimal ordinate, and moving 
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in a clockwise sense, the boundary of a permutomino P can be encoded as 
a word in a four letter alphabet, {N, E, S,W}, where N (resp., E, S, W) 
represents a north (resp. east, south, west) unit step. Any occurence of a 
sequence NE, ES, SW or l^A^ in the word encoding P defines a salient 
point of P, while any occurence of a sequence EN, SE, WS or NW defines 
a reentrant point of P (see Figure [3|). For simplicity of notation and to 
clarify the definition of the construction recalled in Section [3l the reentrant 
points of a convex permutomino are grouped in four classes; in practice, 
the reentrant point determined by a sequence EN (resp. SE, WS, NW) is 
represented with the symbol a (resp. (3, 7, 6). 




NNENNEESEESSWSSWWNWW 

Figure 3: The coding of the boundary of a permutomino, starting from 
A and moving in clockwise sense; its salient (resp. reentrant) points are 
indicated by black (resp. white) squares. 



2.2 Previous enumerative results 

Let us recall the main enumerative results concerning convex permutomi- 
noes. In [SI [9], using bijective techniques, the authors provide enumeration 
of various classes of convex permutominoes, including the parallelogram, the 
directed convex and the stack ones; moreover, a characterization of the per- 
mutations associated with permutominoes of each class is given. Let Cn 
(resp. Vn, T^ni <Sn) be the set of convex (resp. parallelogram, directed 
convex, stack) permutominoes of size n and 

Cn = {MP) ■■ P^Cn}, ^„ = {vri(P) : PeVn}, 
Vn = {7Ti{P) : PeVn}, 5„ = { ^i(P) : P G 5„ } . 
The enumeration results obtained in [HI [9] are shown in Table (H 
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Class of convex 


Number of convex 


Associated set 


permutominoes 


permutominoes with size n 


of permutations 


Parallelogram 


1 'Pn 1 — 


I'Pn 1 — Cn—1 


Directed convex 


W \ — -h 




Stack 


=2"-i 


= 2-1 



Table 1: Enumeration results. 



where c„ is the nth Catalan number, 

1 /2n 

n + 1 \ n 

and 6„ are the central binomial coefficients 




In [7] it was proved, using the ECO method [2], that the number of 
convex permutominoes of size n is: 

2(n + 3)4---(J, n>l. 

The first terms of the sequence are 

1,4, 18,84,394,1836,8468,... 

(sequence A126020 in [H]). The same enumerative result was obtained, 
using a different approach, by Boldi et al. [6] . 

An interesting study of combinatorial properties of the set C„ is in [3]. 

3 ECO construction of convex permutominoes 

In this section we recall the ECO construction of convex permutominoes as 
given in [7]. 

Let Cn be the set of convex permutominoes of size n and let P Cn', the 
number of cells in the rightmost column of P is called the degree of P. Let 
us consider the following properties of a convex permutomino: 
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Ul : the uppermost cell of the rightmost column of P has the maximal 
ordinate among all the cells of the permutoniino; 



U2 : the lowest cell of the rightmost column of P has the minimal ordinate 
among all the cells of the permutomino. 

According to the ECO method [2], it is necessary to define an operator 
■d : Cn ^ 2^^+^ which defines a recursive construction of all the convex 
permutominoes of size (n + 1) in a unique way from the objects of size n. 
The operator {) defined in [7| acts on a convex permutomino performing 
some local expansions on the cells of its rightmost column. Let ci, . . . , c„ 
(resp. ri, . . . , r„) be the columns (resp. rows) of a permutomino P of size n 
numbered from left to right (resp. bottom to top) and let £(ci) (resp. l{ri)) 
be the number of cells in the ith column (resp. row), with 1 < i < n. The 
four operations of i?, denoted by (a), (/3), (7) and {5) are defined as follows: 

(a) if P satisfies condition Ul, then operation (a) adds a new column 
made of ^(c„) + 1 cells on the right of c„, see Figure HI 



(a) 



Figure 4: Operation (a); the added column has been highlighted 



{(3) it can be performed on each cell of Cn] so let di be the ith cell of c^, 
from bottom to top, with 1 < i < i{cn)- Operation (/?) adds a new 
row above the row containing di (of the same length), and add a new 
column on the right of c„, made of i cells, see Figure [5l 

(7) it can be performed on each cell of c„; so let di be the ith cell of Cn, 
from bottom to top, with 1 < i < i{cn)- Operation (7) adds a new 
row below the row containing di (of the same length) , and add a new 
column on the right of Cn made of n — i + 1 cells, see Figure [6j 

(5) if P satisfies condition U2, then operation (6) adds a new column made 
of i{cn) + 1 cells on the right of c„, see Figure [71 
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Figure 5: Operation the cell di is filled in black , the added column and 
row have been highlighted 



(y) 



■3 



Figure 6: Operation (7); the cell di is filled in black , the added column and 
row have been highlighted 



Obviously in any case the obtained permutomino is a convex permu- 
tomino of size n + 1. 

We refer to [7] for further details and proofs. 

4 The generating algorithm 

Our aim is to illustrate an exhaustive generating algorithm for convex per- 
mutominoes, basing on the ECO construction recalled in previous section. 
In the sequel we will refer only to convex permutominoes, simply named 
"permutomino" . 

First of all we define a subset of permutominoes of size n, the so called 
active permutominoes, then we will show the existence of a bijection between 
the active permutominoes of size n and the set of permutominoes of size n— 1. 
Finally, we will define the generating tree of permutominoes of size n. 
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(5) 

> 



Figure 7: Operation (5); the added column has been highhghted 



4.1 Definition of active permutominoes 

A permutomino P is active if the following conditions hold: 

1. the leftmost column contains only one cell, {i{ci) = 1); 

2. the leftmost reentrant point a has abscissa 2. 

That is, the word encoding the boundary of an active permutomino be- 
gins with " NEN . . . and ends with . . . WW" . In Figure [8] are depicted three 
active permutominoes of size 4, while in Figure [9] there are some permutomi- 
noes that, yet having only one cell in the leftmost column, are not active. 



Figure 8: Three active permutominoes of size 4 



Let P be an active permutomino and let g be the row of P containing 
the only cell in the leftmost column (in Figure [5] is highlighted) ; g can 
be the bottom row but it is never the top row. From Proposition [H g ends 
at the same abscissa of either the above and the below row, if it exists; so 
a configuration as that depicted in Figure [10] is not admissible. Therefore, 
if in an active permutomino P of size n we remove g, we obtain a convex 
permutomino P' of size (n — 1). 
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Figure 9: Some not active permutominoes of size 4 



Figure 10: A not admissible configuration 



Let P' be a permutomino of size (n— 1) and let g' the row in P' containing 
the leftmost sahent point with minimal ordinate (point A in Figure [3]). If 
in P' we add, below g' , a row g one cell longer on the left and ending at the 
same abscissa of g' , we obtain an active permutomino P of size re: 

i. if g' is the bottom row, the added row in P' is one cell longer on the left 

and therefore the new permutomino of size n is convex and active (see 
Figure [m a)); 

ii. if g' = Ti, with 1 < z < re, the row rj_i ends at the same abscissa of g'; 

so adding g we obtain a convex and active permutomino of size re (see 
Figure ttUb)). 

This means that it exists a bijection ip between the set Ca^n of active 
permutominoes and the set Cn-i- 

such that ip{P) is the permutomino of size re — 1 obtained by removing g 
from the active permutomino P of size re. 
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a) 



b) 



Figure 11: A permutomino P' of size (n — 1) and the corresponding permu- 
tomino P of size n 



4.2 The exhaustive generating algorithm 

The algorithm we propose for the exhaustive generation of convex permu- 
tominoes of size n is based on the bijection -i/^, defined in Section 14.1^ and 
on the ECO construction of permutominoes recalled in Section [3l 
The generating process is described by an operator (j) so defined: 

(p Operator: 

1. The first permutomino of the generating process is P^, that is the 
permutomino of size n associated with the pair of permutations (tti , 112) 
of [n + 1]: 

vri = (1,2,3, ... ,n,n+ 1) vrs = (2, 3, 4, . . . , n + 1, 1) 
In Figure [l2] is depicted P^ of size 4. 

2. is an active permutomino; let P^ the permutomino of size {n — 1) 
such that P^ = ip{Pn) {P^ is obtained by removing the bottom row of 

n / 

3. Apply operations (7) and (6) of ECO construction to P^ . Every 
application generates a new convex permutomino of size n. 

4. For each new active permutomino Q repeat the following actions until 
active permutominoes are generated: 
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4.1 remove the q row from Q obtaining Q = ip{Q); 

4.2 apply all the possible operations of the ECO construction to Q. 
Every application generates a new convex permutomino of size n. 



Figure 12: -P^ of size 4 



Our strategy can be represented using a rooted tree, say Cn-tree, so de- 
fined: 

1. the root is and it is at level 0; 

2. if <5 € Cn-tree is an active permutomino at level A; > 0, then (j){Q) = 
'^{'^{Q)) is the operator defined in the ECO construction) and every 
P G 4>{Q) is a son of Q and it is at level (/e + 1). For the sake of 
simplicity, we say that P G (l)''+^{P;;). 

In Figure [13] Cs-iree is illustrated. 

Proposition 2 C^-tree contains all and only the convex permutominoes of 
size n, i.e C„-tree = C^. 

Proof. [Only] The permutominoes in Cn-tree are obtained from permu- 
tominoes of size (n — 1) by applying the ECO construction; so, as proved in 
we generate permutominoes of size n. Therefore, since the permutomi- 
noes of size [n — 1) are each other different, the ones of size n are different 
too. 

[All] We must proof that for each permutomino Q (z Cn there exists a path 
from P^ to Q, that is it exists a finite sequence Pq, Pi, . . . , Pk with k € N 
and Pk = Q such that: 

. Pi+i e ^{Pi), o<i<k-i. 
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12 




p ^ 

Figure 14: C2-tree 



In other words, we must prove that there exists k > 1 such that Q € 
(p^{P^). We know that the active permutominoes of size n are as many as 
the permutominoes of size (n — 1). So it is sufficient to proof the following: 

Proposition 3 All the active permutominoes of size n are generated. 

Proof. By induction on the size n. 

Base. For n = 1 there is only the permutomino containing one cell; if 
n = 2 the unique active permutomino is PJ' (see Figure So for n <2 
Proposition [3] yields. 

Inductive hypothesis. Let us assume that all the permutominoes of size 
(n — 1) are generated and let Cn-i-tree be the associated tree. Then, starting 
from the root P^-i it is possible to reach any permutomino of size (n — 1). 
So, for each permutomino Pn-i € Cn-i there exists a k such that: 

Inductive step. Let Pn be an active permutomino of size n and let Qn-i = 
ip{Pn) (Qn-1 is obtained from Pn removing g). So there exists k such that 



Qn-l £ 



But 



so it follows that: 



n-l 



hp: 



that is, there is a path from the root P^ to any permutomino of size (n — 1). 
Therefore, each permutomino of size n is reachable from P^. □ 



13 



4.3 Algorithm cost analysis 

First of all, we will prove that the height of Cn-tree is n. The proof is based 
on the following propositions. 

Proposition 4 Using the ECO construction of Sectionl^ an active permu- 
tomino of size n is generated by one and only one active permutomino of 
size (n — 1). 

Proof. It follows straightforward from the ECO construction which never 
adds a column on the left containing one cell to the permutoninoes. □ 

Proposition 5 Given an active permutomino Q, the longest path starting 
from Q has length j if its a points lie in {2,h),{3,h + l), . . . ,{j + l,h+j — l), 
h being the ordinate of the lefmost a point. 

Proof. The permutomino ip{Q) of size (n — 1) to which the ECO construc- 
tion is applied, is obtained removing from Q the row g, so the leftmost a 
point of Q is removed in ip{Q). Therefore, 4'{Q) will be active, and then, 
from PropositionlU it can generate new active permutominoes, only if it has 
an a point at abscissa 2; thus, Q must have an a point at abscissa 3. In the 
same way, a permutomino generated from ip{Q) will be active only if ■(/'(Q) 
has an a point at abscissa 3, that is if Q has an a point at abscissa 4, and 
so on up to {j + 1). □ 

The permutomino has {n — 1) consecutive a points in (2, 2), ... , (n, n), 
so, from Proposition El the longest path starting from has length {n — 1). 
Therefore, the height of Cn-tree is n. 

Prom the generating alghorithm of Section 14.21 it follows that the gen- 
eration of a permutomino P at level k in the Cn-tree depends only on a 
permutomino at level {k — 1), that is, Cn-tree is generated level by level. 
Therefore, since the generation of a permutomino of size n from one of size 
— has a constant cost, we may conclude that the cost of the exhaustive 
generating algorithm is proportional to the number of permutominoes of size 
n. 
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